Exchange interaction and Fano resonances in diatomic molecular systems 
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We propose a mechanism to use scanning tunneling microscopy (STM) for direct measurements 
of the two-electron singlet-triplet exchange splitting J in diatomic molecular systems, unsing the 
coupling between the molecule and the substrate electrons. The different pathways for electrons 
lead to interference effects and generate kinks in the differential conductance at the energies for 
the singlet and triplet states. These features are related to Fano resonance due to the branched 
electron wave functions. The ratio between the tunneling amplitudes through the two atoms can be 
modulated by spatial movements of the tip along the surface. 

PACS numbers: 73.23.-b, 68.37.Ef, 72.15.Qm 



There are various techniques that allows one to detect 
and manipulate spin states in the solid state, which at- 
tract lots of interest. A partial list include optical detec- 
tion of electron spin resonance (ESR) in a single molecule 
[H, tunneling through a quantum dot [§], and, more re- 
cently, ESR-scanning tunneling microscopy (ESR-STM) 
technique 0, [|]. The interest in ESR-STM is due to 
the possibility of manipulating single spins H, H, 0], 
something which is crucial in spintronics and quantum 
information. Experimentally, modulation in the tun- 
neling current has been observed by STM using spin- 
unpolarized electron beam 0, Q . Lately, there has also 
been a growing interest in using spin-polarized electron 
beam for direct detection of spin structures [8|], as well 
as utilizing the inelastic electron scanning tunneling spec- 
troscopy (IETS) for detection of local spatial variations 
in electron-boson coupling in molecular systems [13] • 

Typically in STM measurements with an object lo- 
cated on a substrate surface, the tunneling current can 
either go directly between the STM tip and the substrate 
or go via the object. The tunneling electrons are thus 
branched between different pathways, which gives rise to 
interference effects when thepartial waves merge into one 
in the tip or the substrate [9(. This interference leads to 
a suppressed transmission probability for the tunneling 
electrons at certain energies. The suppressed transmis- 
sion is a fingerprint of Fano resonances , and generally 
appear in systems where tunneling electrons are branched 
between different pathways. Recently, Fano resonances 
have been studied in double and triple quantum dots sys- 
tems 12j, where the different pathways are constituted 
of the different quantum dots. 

In this Letter we propose a method to measure the 
two-electron singlet-triplet (S-T) exchange splitting J in 
a diatomic molecule by means of STM. The presence of 
two pathways for the tunneling current between the tip 
and the substrate, through the diatomic molecule, gives 
rise to interference effects (Fano resonance) between the 
electron waves traveling through the singlet and triplet 



states. In the direct tunneling between the tip and sub- 
strate via the molecule, the probability for the tunneling 
is proportional to To- In addition, because of the phase 
space branching of tunneling possibilities, the tunneling 
probability is multiplied with the interference probabili- 
tiy 7 2 , hence, the characteristic energy width of the anti- 
resonance is 7 2 ro. Clearly, the anti-resonances will be 
measurable in the second derivative of the current when- 
ever 7 -c 1. 

Fano resonances can be realized in a variety of sys- 
tem, ranging from systems with interactions between 
continuum states and a localized state, to systems where 
the branching of the wave function through diatomic 
molecules. In the case we consider here, we also have 
to include the fact that the one-electron states in a two- 
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FIG. 1: (Color online) a) Cartoon of the diatomic A + B 
molecule coupled to the tip and substrate. The particles 
interact via exchange interaction J. b) Transmission for 
transitions between the singlet/triplet and the anti-bonding 
(dashed) and bonding (solid) one-electron states, c) Phase 
space branching of the tunneling wave functions between the 
singlet \S) and triplet \T) states. The respective amplitudes 
of the tunneling through these states have to be added, re- 
sulting in the Fano like features. 
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level system with the levels being resonant, consist of 
an anti-bonding and a bonding state (both being spin- 
degenerate). This modifies the expected transport prop- 
erties such that transitions between the triplet and the 
anti-bonding one-electron state generate a dip in the 
transmission at the energy for this transition, see Fig. 
[T]b) (dashed). On the other hand, transitions between 
the singlet and the bonding one-electron state give rise 
to a dip in the transmission at the energy for this tran- 
sition, see Fig. [1] b) (solid). In turn, the features that 
appear at voltages corresponding to both the singlet and 
triplet states, provide a unique fingerprint which allows 
the read out of the S-T exchange splitting J. 

In a regular Fano resonance, interference occur be- 
tween the different tunneling paths in real space, one path 
going through the molecule to the substrate whereas the 
other goes directly into the substrate. Here, on the other 
hand, the interference occur between different pathways 
in phase space where the tunneling can occur through 
a singlet \S) or triplet \T). The respective amplitudes 
of the tunneling through theses states have to be added, 
resulting in the Fano like features, see Fig. Q]c). 

In general, Fano resonance appearing at voltages that 
correspond to both the singlet and the triplet states is 
expected to be found in any model which includes one- 
electron states and two-electron singlet and triplet con- 
figurations, and which accounts for the S-T exchange 
splitting J. For simplicity, however, we will perform our 
calculations for the interference between the singlet \S) 
and triplet \T) = \S = 1, S z = 0). We checked that 
the interference arising from the other triplet state \S = 
1, S z — ±1) leads to results which only renormalizes the 
coefficients in the final expressions for the transmission, 
and therefore can be omitted. Hence, we use a simplified 
model for the molecule, which both provides the bonding 
and anti-bonding one-electron states, as well as the two- 
electron singlet \S) and triplet \T) = \S = 1, S z = 0) 
states, i.e. 



d Ba d Ba ) + J{S+S B +S 



A B 



(1) 

Here, ^ao/Bg (^Aa/Ba) creates (annihilates) an elec- 
tron in atom A/B at the energy sq, whereas J de- 
notes the exchange splitting, and S A = d A +d,Ai and 

S A = d A ,d,A-[, and analogously for B. In this model, 
the one-electron states are the (spin-degenerate) anti- 
bonding and bonding states | lcr) = (d Aa — d' Ba )\0) / y/2 
and 1 2a) = (d Aa + d' Ba )\0) /y/2, respectively, whereas the 
two-electron states are \S/T) = [d B >d A ^d B +d A ,]\Q)/"\/2 
for the singlet (S) and triplet (T), respectively. In order 
to keep the number of parameters to a minimum, we ig- 
nore the splitting between the bonding and anti-bonding 
one-electron states. This assumption is not crucial for our 
subsequent analysis, since the Fano resonance we discuss 



arise due to interference between the singlet and triplet 
states, and not between the one-electron states. In order 
to make contact with experiments, we assume that J cor- 
responds to the exchange splitting which is renormalized 
by the coupling to the substrate, that is, the J that will 
be detected by any measurement in a given setup. 

To illustrate the electron interference along different 
paths consider an electron tunneling from the tip to 
the substrate via the molecule. Because of the different 
paths, the electron acquires a different phase depending 
on whether it travels through A or B and, in addition, 
whether it makes a transition between the triplet or sin- 
glet state and the bonding or anti-bonding one-electron 
states. To be specific, suppose a spin f electron leaves 
the molecule from the singlet state to a medium that 
hybridizes equally strong to A and B, e.g. 

{d A] + d BT )|5) = (-4 4 - d Al )\0)/V2 = -|2 |). (2) 

The final state is orthogonal to the anti-bonding state 
|1 J.) Hence, the singlet state couples to the bonding one- 
electron state. Likewise, supposing that a spin f electron 
tunnels out from the triplet state, e.g. 

(d A] + d B] )\T) = (-d Bi + d Al )\Q)/V2 = \li), (3) 

shows that the triplet only couples to the anti-bonding 
one-electron state. 

The tunneling between the tip/substrate and the 
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VBkaC ka dBa- + H.c], where d_ creates an electron in the 



molecule is modeled by Ht 

iT y-J „ I f-f rt 1 titVi at-a 

tip/substrate, whereas v A k a and VBka are the tunneling 
rates between the tip/substrate and A and B, respec- 
tively Reformulating the molecule in terms of its eigen- 
states, e.g. H AB = E CT ,„=i,2 E n \na)(na\ + E S \S)(S\ + 
Et\T)(T\ with the eigenenergies E n — Eq, n — 1,2, cor- 
responding to the states \na), and E S / T = 2e T «//2, 
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FIG. 2: Example of tunneling events leading to the inter- 
ference. First a) an electron leaves the molecule through a 
transition to the anti-bonding one-electron state \2a). Sec- 
ond, an electron enters the molecule through a transition to 
either the singlet state b) or to the triplet state b'). 
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and neglecting the transitions between the empty and 
one-electron states, the tunneling Hamiltonian becomes 



whereas 



lX [ v k<mxc , k(r \na 



}(X\ + H.c], where X = S,T, 



VkanX = v Akr7 (na\d A a\X) + v Bka {7ia\d B a\X) 



(4) 



This form of the tunneling rate reflects the real space 
branching of the electrons tunneling between the tip and 
the substrate via the molecule. From this description of 

L I R 

the tunneling rate we define the couplings r L, , = 
2nJ2keL/R v k<TnX v k<TriX'8(u ~ Ska), where L/R denote 
states in the tip/substrate. In the present study, the spin 
currents through the system are identical, thus we do not 
explicitly indicate the spin in the couplings. 

Suppose that the molecule is in the two-electron sin- 
glet state. The tunneling current between the tip and 
the substrate is then mediated by sequences like \S) — > 
\na) — > \X), where X either of the singlet (S) or the 
triplet (T) states. Hence, there occur indirect S-T transi- 
tions, in the sense that the tunneling may cause a singlet 
configuration to be turned into a triplet, and vice versa, 
see Fig. [3] These indirect S-T transitions give rise to 
the interference effects between the different transport 
channels, which implies that multiple scattering events 
have to be taken into account in the description of the 
tunneling through the molecule. These tunneling events 
leading to the interference resemble cotunneling often dis- 
cussed in this context Here, the multiple scatter- 
ing events are accounted for through the Green functions 
(GFs) G„xx>n{t, 0) [14], where subscripts nX denote the 
transitions \na)(X\, whereas subscripts Xn denote the 
conjugate transitions |X)(na|. 

We assume, for simplicity, that both atoms couple 
equally to the substrate and model this by letting VAqa — 
VBqa = Vqa for qa € R. From Eqs. ([3]) and ((H this 
leads to a vanishing probability for transitions between 
the triplet and anti-bonding one-electron state, since 
|*V2t| = \v qa \\{(2a\d Aa \T) + (2a\d Ba \T})\ = \v qa \\ - 
1/2 + 1/2 1 = 0. Likewise, from Eqs. © and (gj), we 
note that the probability for transitions between the sin- 
glet and bonding one-electron state vanishes. Hence, the 
only non- vanishing coupling matrix elements between the 
molecule and the substrate are given by r^ T = Tt/ S = 
To/2, where T = ^J2k^L,R Wa\ 2 5{uo - e fc(J ). As we 
mentioned, the Fano resonance feature should be present 
for any bonding and anti-bonding splitting, since the in- 
terference we discuss occur between the different tunnel- 
ing processes \S) — > \na) — > \S) and \S) — ► \na) — > |T), 
and not between the one-electron states, see Fig. [51 

The tunneling between the tip and the molecule, how- 
ever, is assumed to depend on the spatial position of 
the tip relatively A and B. Introducing the functions 
1a/b{\ y ~ v a/b\), we model the spatial dependence by 
putting VApa = 7A(|r - ta\)v P(X and VBpa = 7s(|r - 
r.B|)i>p<T for pa G L. The coupling matrix between 



the anti-bonding one-electron state and the two-electron 
states for the tunneling between the tip and the molecule 
then becomes 



rf 



1 1T 



1 1TS1 



1 1ST1 1 IS 

{lA+lB? 1a 
1 1 (1A-1B? 



7l 



T A -7b 



(5) 



The coupling matrix between the bonding state and the 
two-electron states is obtained by letting ± — ► =p in the 
diagonal elements in Eq. ([5]). The assumption that the 
atoms couple equally to the substrate is not crucial. Un- 
equal couplings between the atoms the substrate would 
result in the coupling matrices Ty 2 being full (struc- 
turally equal to Ty 2 ). However, since the tip is mov- 
able, inequalities in the coupling between the atoms and 
the substrate can be canceled by relocating the tip to a 
position where it couples equally to both atoms. 

Using the equation of motion technique we find that 
the Fourier transformed retarded GF, to a good approx- 
imation, can be solved by flBj 



C„(w) 



C„(w) ' 



3£sM = E 



\Vk<rnS\ 



-'nTSn 



tO - Ska + M + ' 

V kanT V k<JnS 
LU ~ Ekcr + M + ' 



, (6a) 
(6b) 
(6c) 

(6d) 



where E s - E n = e - J/2. The GFs G r nSSn and G r nSTn 
are obtained from Eq. © by letting T (S) -> S (T), 
where Et — E n = e + J/2. In Eq. © we have used the 
function C n (uj) = det G^ -1 = (u — u) n +)(oj — w n _), that 
is the poles of the GF, which are given by 



(7a) 




J- «^(1 + 1a7b) 



. 2± = £0 -^( 1+ ^i 



r 7i-7^ 2 



(7b) 



±- 



J + i — (1 + 1a1b) 



r 7i - Yb 



2 \ 2 



These expressions are obtained by neglecting the 
real parts of the self-energies S^ TSn , we have 



J l(2 



-i{T[ 



(2) 



rf (2) )/2. Here, subscript 1 (2) 



refers to the anti-bonding (bonding) one-electron state, 
whereas + (— ) refers to the triplet (singlet) two-electron 
state. Hence, the energy loi + corresponds to transitions 
between the triplet and anti-bonding states. 
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From Eqs. ([7]) we extract the widths of the transi- 
tions between the one- and two-electron states. Assum- 
ing 7s = k7a, where ja "C 1 and < k < 1, we find 
that the widths of the transitions between the triplet (sin- 
glet) and bonding (anti-bonding) states are significantly 
smaller than the transitions between the triplet (singlet) 
and the anti-bonding (bonding) states. Indeed, while the 
widths of the main peaks at uj\ + and 102- are approxi- 
mately To/4, the widths of the narrow peaks at wi_ and 
102+ are approximately (r /4)(l + k 2 )7^/4. Hence, for 
small we expect sharp features in the differential con- 
ductance for the system. However, due of the small width 
of the dips, these features will be easier to observe in the 
derivative of the differential conductance, i.e. the second 
derivative of the current / ~ tr J T [ui)[f l{uS) — fn(uj)]duj 
which here is given by 



d 2 i 
W 2 



/ tanh [I3{lo - Ml )/2] 
Vcosh 2 \ft(u) - fi L )/2] 

tanh \fl(u - Mr)/2] 



cosh 



dw, (8) 



where T n {u>) 
sion coefficient 



tr r^G^(w)r^G^(w) is the transmis- 
[ljj , whereas n L / R is the chemical poten- 
tial of the tip/substrate. Under the given conditions we 
find that the total transmission T = ^ n 7~ n through the 
molecule is given by 



-j-) (7-4 + 1b) 2 



LU — 



gp + J/2 



LU — 



£0 - J/2 



C 2 (w) 



1,(9) 



which shows dips at the energies corresponding to the 
transitions between the one-electron and two-electron 
singlet and triplet states, e.g. at lo — e T J/2, respec- 
tively, independent of Ja/b- From the above analysis 
of the widths, these dips are expected to be sharp since 
their positions in the transmission exactly correspond to 
the roots lo\- and L02+, respectively, for small 7^4. 

Asymmetric coupling between both the tip and sub- 
strate and the atoms generate a shift in the positions of 
the transmission dips [17J . This result shows that the dis- 
tance between the transmission dips deviates from J by 
at most 6 % for asymmetries VA P ( q )cr/vBp{q)cr = 1l{r) < 



0.7 17J. Thus, even for asymmetric coupling, this will 



still be a good measure of J. 

In conclusion we propose to use the interference be- 
tween different paths of the tunneling electron wavefunc- 
tion as a tool to detect the S-T splitting in a single 
molecule. This splitting could be observed as features 
in the tunneling conductance. Becasue the features we 
find are rather narrow they can be, more prominently, 
revealed in the d 2 I/dV 2 peaks that would correspond 
to the S-T splitting. The splitting produces interference 
between the possible pathways for electrons to tunnel 



through virtual states with different ener gies 
to the classical Fano resonance arguments 11 



identical 
This ef- 
fect is similar to IETS features seen in inelastic scattering 
off non-magnetic and magnetic excitations in molecules 
[lfj| . With current expertimental capabilties in STM one 
can easily address single diatomic molecules with either 
spin-polarized STM [8j or with position dependent IETS 
measurements 10]. 
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